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Abstract 

Migration is a fundamental stage in the life history of several taxa, including birds, and is under strong selective pressure. At 
present, the only data that may allow for both an assessment of patterns of bird migration and for retrospective analyses of 
changes in migration timing are the databases of ring recoveries. We used ring recoveries of the Barn Swallow Hirundo 
rustica collected from 1908-2008 in Europe to model the calendar date at which a given proportion of birds is expected to 
have reached a given geographical area ('progression of migration') and to investigate the change in timing of migration 
over the same areas between three time periods (1908-1969, 1970-1990, 1991-2008). The analyses were conducted using 
binomial conditional autoregressive (CAR) mixed models. We first concentrated on data from the British Isles and then 
expanded the models to western Europe and north Africa. We produced maps of the progression of migration that 
disclosed local patterns of migration consistent with those obtained from the analyses of the movements of ringed 
individuals. Timing of migration estimated from our model is consistent with data on migration phenology of the Barn 
Swallow available in the literature, but in some cases it is later than that estimated by data collected at ringing stations, 
which, however, may not be representative of migration phenology over large geographical areas. The comparison of 
median migration date estimated over the same geographical area among time periods showed no significant 
advancement of spring migration over the whole of Europe, but a significant advancement of autumn migration in 
southern Europe. Our modelling approach can be generalized to any records of ringing date and locality of individuals 
including those which have not been recovered subsequently, as well as to geo-referenced databases of sightings of 
migratory individuals. 
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Introduction 

Migration is widespread in nature and several taxa, from insects 
to fishes, amphibians, birds and mammals, undertake annual 
"incredible journeys" that represent key stages in their yearly cycle 
[1—3]. Being able to fly, birds are the taxon where migratoriness is 
most widespread, and on which the majority of migration studies 
have focused [2]. As a fundamental feature in the life-history of 
birds, migration is under strong selective pressures [1]. However, a 
large amount of genetic variability in migratoriness, timing of 
migration, and migration strategies exists in bird populations. In 
addition, individuals show a high degree of phenotypic plasticity in 
migration strategy [4]. Genetic variability and phenotypic 
plasticity allow birds to adjust their migration strategies according 
to changes in climate and ecological conditions. Indeed, changes 
in the timing (phenology) of migration are considered signals of the 
impact of current climate changes on the biosphere [5-7]. 



However, many species of migrant birds are declining, probably 
because they are not able to sufficiently adjust the timing of their 
annual life-cycle to match new climatic conditions [8-10]. 

Studies of bird migration are hampered by the difficulty of 
tracking small-sized species, which represent the large majority of 
migratory birds. New miniaturized and cheap technological 
devices, like light-level geo-locators, are bridging this gap in our 
knowledge [1 1] although the information they provide must be 
interpreted with caution because of the impact they may have on 
individual fitness [12-14], and, potentially, also on migration 
timing and routes [12]. Most importantly, these data, together 
with stable isotope analyses of museum specimens [15], will likely 
represent, for a long time, the only sources of information allowing 
for retrospective analyses of changes in bird migration strategies 
through time. 

The analysis of ring recoveries is hampered by several 
difficulties, the main ones being the large spatial and temporal 
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heterogeneity in ringing effort and in the probability of recovery of 
a ringed individual [16]. Nevertheless, these data have been useful 
to study large scale patterns of individual distribution, like 
migratory connectivity, and long term variation in bird distribu- 
tion [17—20]. Retrospective analyses of ring recoveries may also 
allow us to quantify both the progression and the timing of 
migration. Progression of migration is defined here as the 
proportion of individuals of a given migratory population that, 
at any given time, have reached or have passed over a given place 
during migration, while migration timing is defined here as the 
date at which a given proportion of individuals have reached a 
given location. 

In the present paper we show how ring recoveries collected 
throughout Europe and North Africa in 1908-2008 and stored in 
the EURING databank (EDB, www.euring.org) can be used to 
model the progression and the timing of migration of bird species, 
and its variation over time, using a small passerine, the Barn 
Swallow Hirundo rustica, as an example. 

Progression and timing of migration can be statistically 
modelled by fitting the complementary log-log ('cloglog' hereafter) 
function [21], which is very similar to the logistic function (see 
Text SI), to the cumulated proportion of individuals that have 
reached a given site by a given date. This modelling approach has 
the advantage that parameters of the interpolated cloglog curve 
describe both the progression and the timing of migration at a 
given site, since the function allows us to estimate both the 
expected proportion of migrants that have arrived or passed at a 
given date, and, by model inversion, the date when a given 
proportion of migrants is expected to have arrived at or passed 
over a given site. 

We first modelled progression of spring migration over the 
British Isles, taking advantage of the very large amount of data 
available for this area. We then tentatively extended the same 
model to Europe and north Africa, where data are sparser. Paucity 
of data from eastern Europe and the Middle East forced us to 
restrict the analysis to western Europe and the western part of 
north Africa (western Europe and north Africa hereafter). The 
same approach was also used to model autumn migration, first in 
the data-rich British Isles and then in western Europe and north 
Africa. Secondly, we modelled the variation of migration 
phenology over time. To this aim we divided the dataset into 
three periods containing approximately a similar amount of data 
(1908-1969, 1970-1990, and 1991-2008) and compared median 
migration date estimated by the cloglog functions describing 
progression and phenology of bird migration in different 
geographical areas. Also in this case analyses were run separately 
for the British Isles and for western Europe and north Africa and 
for spring and autumn migration, respectively. 

Materials and Methods 

Datasets 

For individually ringed birds, the EDB includes information on 
date and locality at ringing, as well as at any subsequent 
encounter. These data will hereafter be defined as "ring 
recoveries" in order to distinguish them from records of ringing 
date and locality of individuals, which have not been subsequendy 
recovered ("ringing data" hereafter). Hence ring recoveries 
include both ringing and finding information of any bird that 
has been re-encountered. 

The EDB almost exclusively include ring recoveries (a few 
ringing data has been recently included in the EDB, but they were 
not considered in the present study). Before the analyses we 



carefully checked the consistency of data in our dataset and 
excluded any dubious data (details not shown). 

The datasets used to model migration in the British Isles consist 
of 1983 ring recoveries during spring migration (March-June) and 
8429 ring recoveries during autumn migration (August-October), 
while those used for Western Europe and North Africa consist of 
11918 ring recoveries (including ring recoveries from the British 
Isles) collected during spring migration (February-June) and 28832 
during autumn migration (August-November) (figure SI). Periods 
of spring and autumn migration were chosen according to Cramp 
[22]. Since the relevant information for this analysis is the date at 
which an individual was observed in a given geographical location, 
we used all records of individuals found either alive or dead, and 
retained repeated records of the same individual. Records of 
nestlings i.e. birds ringed at the nest and unable to fly (EURING 
age code equal to 1; see Speek et al. [23] for further details on the 
EURING exchange code 2000) were excluded, as well as records 
of individuals found dead, but not fresh (EURING code 
'condition' either equal to 1 or 3 [23]) and those of individuals 
whose recovery date is known with an accuracy larger than 3 days 
either side of the reported date (EURING code 'accuracy of date' 
in 0-2 [23]). Date of recovery is given in days with January 1 as 
day 1. 

Conditional autoregressive models 

The British Isles were divided in 38 cells of 1 .5° x 1 .5° latitude x 
longitude, while western Europe and north Africa were divided 
into 67 4°x4° cells, and each recovery was assigned to a cell 
(Figure SI). Only cells with ring recoveries recorded in at least four 
different dates were included in the analyses because interpolation 
of cloglog curves requires at least four data points per cell. Cell size 
was chosen as to maximize the geographical coverage of cells 
suitable for analyses. Only cells to the north of latitude 26° N and 
to the west of longitude 26° E were considered because data were 
too scattered outside this area. Some cells could not be included in 
all analyses due to paucity of data. For the British Isles, 27 cells 
were included in the spring analyses and 29 in autumn, while for 
western Europe and north Africa there were 59 cells in spring and 
53 in autumn. 

The analytical procedure interpolated the cumulated propor- 
tion of Barn Swallows recorded in a cell at each date over the 
periods of spring or autumn migration, by also accounting for the 
spatial autocorrelation of data recorded at the same time in 
adjacent cells. Let nj t be the cumulated number of individuals 
observed in cell j until date t, irrespective of the year of recovery, 
T 

Nj = nj t be the total number of Barn Swallows in cell j and 
t=i 

Pj t = rij t /Nj the proportion of Barn Swallows recovered until date 
t, t = 1 , . . . , T, T being the end of the period of interest (Figure S2). 
All recoveries were used irrespective of year as data were sparse for 
some cells. Ordinary binomial regression can be adopted to 
estimate the cumulative proportion of arrivals in cell j at any given 
date as a function of a set of secondary variables. We modelled the 
occurrences in a cell as a linear function of the date on a cloglog 
scale since this scale is the most appropriate to model spatial point 
patterns on a geographical grid [21]. To account for spatial 
autocorrelation and avoid biased estimates, we specified an 
autobinomial spatial model for arrivals by including among the 
linear predictors a spatial covariate obtained by calculating the 
weighted average proportion of Barn Swallows that, at any given 
date, had reached the cells immediately adjacent to any given cell. 
To account for potentially different cell counts, each cell in the 
neighbourhood was weighted by the proportion Wk = N k /Naj of 



PLOS ONE | www.plosone.org 



2 



July 2014 | Volume 9 | Issue 7 | e102440 



Modelling Migration by Ring Recoveries 



the arrivals at any cell k in the neighbourhood 3j of cell j out of the 
global number of arrivals in the neighbourhood, i.e. Ngj = Nk. 

Cells were considered adjacent when they shared a side or a 
vertex ('queen' configuration [24]). To account for the inter-cell 
variability in patterns of migration through time, cell identity was 
entered as a random grouping factor and date as a random slope 
at the cell level. More formally the model is specified by 

log ( - log ( 1 - E(p Jt )) ) = a + (3t + y ^ w kPkt + Aj + Gj t 

where Aj~jV(0, T 2 ) and Gj~N(0, u 2 ), obtaining a Conditional 
Autoregressive (CAR) binomial Generalized Linear Mixed Model 
(GLMM). 

Note that this model accounts for different numbers of 
observations at each cell in two ways. Firstly, the dependent 
variable pj t is specified in the model as the ratio between the 
cumulative number of individuals that have reached cell j until 
time t over the total number of individuals in the cell (actually n it / 
(Nj - n it ) in the procedure we used for analyses [25]). The variance 
of the dependent variable is therefore calculated by taking into 
account the total number of observations at a cell, thus giving 
larger weight to cells with more observations. Secondly, the spatial 
autocovariate is also calculated by giving larger weight to cells in 
the neighbourhood with more observations. 

The ability of the model to correctly interpolate the observed 
proportion of Barn Swallows that had arrived or had migrated 
over a given cell in a given date was estimated by calculating a 
pseudo R 2 equal to the squared correlation coefficient between 
observed and estimated cumulated proportions at each cell and 
date for which there were observations (Efron's pseudo R 2 for 
binomial models; [26]). This is a measure of the predictive ability 
of the model, similar to common R 2 of linear models, which is 
undefined for binomial models [27]. 

These analyses were performed by the Imer procedure in the 
lme4 package [25] in R 2.15.2 [28]. The cloglog function was 
interpolated by specifying the cloglog link function in the Imer 
procedure. 

Model inversion and map production 

The CAR model was used to predict the exact date at which a 
given proportion of Barn Swallows had been recovered in a given 
cell and used to produce contour maps of the date at which a given 
proportion of Barn Swallows have reached a particular geograph- 
ical location. In particular, the date t at which a proportion p of 
swallows are estimated to have reached or have migrated over a 
given cell j can be calculated as 

1 j = (log(-log(l-p))-«j)/p i 

where otj and fSj are, respectively, the values of intercept and slope 
estimated by the CAR binomial GLMM at cell j. 

In this paper we produced maps of the calendar date at which 
15%, 50% and 85% of Barn Swallows have reached a given cell 
during spring migration, and maps of the date at which the same 
percentage of Barn Swallows was still in the cell during the period 
of autumn migration. In addition, we estimated arrival date of the 
5% of Barn Swallows, and date when 5% of Barn Swallows were 
still in the cell for comparison with arrival dates of the first 
individual and departure date of the last individual observed at 
ringing stations or other locations. 



Contour maps allow speculations on migration flyways as the 
contours are isochrones that connect geographical areas showing 
the same phenology. If Barn Swallows follow flyways during their 
migration, geographical localities along the flyways might be 
reached by a given proportion of Barn Swallows earlier in the 
season than the surrounding areas where the migration movement 
is less intense. Hence, we expect map contours to show a reverse-U 
shape in the flyway direction. 

Consistency of observed and model-predicted 
phenology 

We aimed at comparing our model-predicted estimates of 
migration phenology with known information of phenology 
derived from the literature. We considered both quantitative 
estimates of first and mean/median arrival dates or departure 
dates of the last individual from time series of ringing/ observation 
and qualitative descriptions of migration phenology (Table SI). 
Quantitative phenological data from time series at a given 
geographical location were compared to arrival/ departure dates 
estimated by our model for the corresponding percentage and cell. 

Qualitative descriptions of migration phenology were also 
entered in the analysis by converting them to a quantitative 
estimate (Table SI). We acknowledge that this procedure is based 
on a subjective interpretation of the qualitative description, but 
note that excluding these data from the analyses did not alter the 
results of the following analyses (details not shown). 

Consistency and agreement of phenological estimates from our 
models and observed phenology were assessed by calculating the 
repeatability [29;30] between phenological estimates from the 
literature and those estimated by our models. Due to paucity of 
data, we pooled data from both autumn and spring migration. To 
avoid unduly inflating repeatability due to the (obvious) difference 
in dates among spring and autumn migration we used the 
following procedure to centre the data before the analysis. We 
calculated the mean value of both observed and estimated values 
for spring migration (common mean for spring migration), and the 
mean value of both observed and estimated values for autumn 
migration (common mean from autumn migration). We then 
subtracted the common mean for spring migration from both 
observed and estimated values for spring migration, and the 
common mean for autumn migration from both observed and 
estimated values for autumn migration. 

Maps of ring recoveries 

Maps of ring recoveries assist with the interpretation of 
movement patterns inferred from contour maps by showing the 
actual movement of individuals. Information on the movements of 
individuals is included in ring recoveries, but was not used in the 
analysis on which contour maps are based. The only information 
necessary to produce contour maps is indeed the date and the 
position where an individual has been observed. 

Maps of ring recoveries were produced by connecting the 
positions where an individual was observed, irrespective of the 
year of recovery. Only records during spring or autumn migration 
were included. As maps of ring recoveries were used for 
comparison with contour maps, they only included individuals 
that moved a range of distances comparable to those that could be 
inferred from contour maps. For this reason, only individuals that 
moved between 1 and 8 degrees of latitude or longitude were 
included in maps of ring recoveries. Indeed, lines connecting the 
positions of individuals recovered at longer or shorter distances 
only complicate these maps without providing useful information 
on patterns of migration. 



PLOS ONE | www.plosone.org 



3 



July 2014 | Volume 9 | Issue 7 | e102440 



Modelling Migration by Ring Recoveries 



Temporal variation in migration phenology 

We investigated whether median migration dates varied over 
time and whether changes in the timing of migration differed 
among geographical areas [31]. To this end, the dataset was 
divided into three time periods, 1908-1969, 1970-1990, and 
1991-2008. Time limits for these periods were chosen to include a 
similar amount of data for each period (details not shown). In 
addition, the British Isles were divided in two latitudinal belts 
north or south of latitude 53° 45' N (Figure SI), while western 
Europe and north Africa were divided into three latitudinal belts 
(northern Europe: >50° N, central Europe: >42° N & <50° N, 
southern Europe and north Africa: <42° N; see Figure SI). These 
thresholds were chosen so as to have at each belt a sufficient 
number of cells for statistical analyses in each period. For each cell 
the cumulative proportion of Barn Swallows that had been 
recovered by a given date within each period was then modelled as 
a cloglog function of date. The date when the median Barn 
Swallow was expected to arrive at any cell ('median migration 
date' or to hereafter) in each period was calculated from the fitted 
cloglog curve (t 0 = (log( — log(0.5)) — oc)/P where a and P are, 
respectively, the values of intercept and slope of the (linearized) 
cloglog curve fitted at each cell). 

Median migration dates in each period and cell were then 
analysed by using CAR linear mixed models whereby period 
(three-level factor), belt (two-level factor in the analyses on the 
British Isles and three-level factor in those on western Europe and 
north Africa) and their interaction were entered as predictors 
together with a spatial autocovariate, while cell was entered as a 
random grouping factor. The value of the spatial autocovariate 
Wj,t for cell j in each period n was here calculated by averaging 
median migration dates \x^ n at period n of the cells in the 
neighbourhood 3j„ of cell j. Formally: 



We also corrected our models for heteroscedasticity because 
graphical exploration showed that variance in median migration 
dates differed widely between periods. In these models, a 
significant belt by period interaction would indicate that changes 
in timing of bird migration between periods differed according to 
the geographical position of cells. 

These analyses were performed with the Ime procedure in the 
nlme package [32] in R 2.15.2. Models were corrected for 
heteroscedasticity by specifying an among-period varldenl weight- 
ing function in the Ime procedure [32]. 

Results 

Maps obtained from CAR models 

British Isles. Maps of spring migration show a general 
northward progression of migration, but also show some local 
patterns (Figure 1 A). The maps of spring migration in the British 
Isles indicate that Barn Swallows arrive earlier in central Ireland 
(estimated migrtion date of the first 15% of Barn Swallows at cell 
El is 109 — 19 April; see Figure SI A for cell IDs) than in the rest of 
the British Isles (data from southern and western Ireland were 
unavailable). They then appear to move northwards toward south- 
west Scotland, where the first 15% of Barn Swallows is expected to 
pass around 8 May (date = 128). 

South-west England is reached by the first 15% of Barn 
Swallows at the beginning of May (124 = 4 May), then Barn 



Swallows seem to move northwards in two main directions, on the 
one side toward Wales and north-west England, and on the other 
along the western coast of northern England. Barn Swallows arrive 
latest in northern Scotland and the Orkney Islands (15% 
approximately at 140 = 20 May). The bulk of migration (50%) 
transits in Britain around 30 May (150), with the only exception of 
central Ireland, where it is earlier (128 = 8 May), and of northern 
Scotland, where it is later (154= 3 June). The last Barn Swallows 
(85%) pass through central Ireland on around 30 May (85% is 
estimated at date =150 at cell El), then probably move north 
towards Scotland, which is reached by the last 85% of Barn 
Swallows on around 13 June (164; Figure 1A). 

The CAR mixed model used to produce the map of spring 
migration over the British Isles interpolated the observed 
proportion of Barn Swallows at each cell in each date with great 
accuracy (R 2 = 0.95). In addition, the pattern of migration 
depicted above is consistent with the movements of individual 
Barn Swallows documented by ring recoveries. Indeed, Figure 2A 
shows movements of individuals between the Channel Islands and 
Ireland, thus suggesting direct movements of Barn Swallows 
towards Ireland. 

Maps of autumn migration in the British Isles showed a reverse 
pattern, with Barn Swallows moving south-east from Wales and 
Scotland. In addition, Barn Swallows seem to move through the 
western part of Ireland earlier than from Scotland, and leave 
eastern England last (Figure IB). Migratory movements are more 
synchronous during autumn than spring migration, as indicated by 
the lower maximum difference in dates represented by isochrones 
on maps of the autumn compared to spring migration. In addition, 
autumn migration routes seem to follow a more eastward 
direction, as suggested by the shape of isochrones which point 
east in southern England. The predominant eastward movement 
during autumn migration is confirmed also by the maps of ring 
recoveries (Figure 2B). More detailed patterns are difficult to assess 
in this map, probably due to the synchrony of movements during 
autumn migration. The model used to produce the map 
interpolated the observed proportion of Barn Swallows at each 
cell in each calendar date with great accuracy (R 2 = 0.97). 

Western Europe and North Africa. Maps of spring 
migration over western Europe and north Africa (Figure 3A) were 
based on a model that interpolated the observed data with great 
accuracy (R 2 = 0.96). They show an early transit of birds during 
spring migration in the Iberian peninsula, with the first 15% of 
Barn Swallows in southern Portugal on 1 March (60) and in 
central Spain and southern France on 10 April (100). They seem 
then to spread north-eastwards in France and the rest of northern 
Europe, reaching southern Sweden on 20 May (140). 

The first 15% of Barn Swallows reaches southern Italy at the 
beginning of April, then they move towards the Balkans. 
Interestingly, the maps of the 15%, 50% and 85% of Barn 
Swallows suggest a progressive eastward shift of the northward 
turn of isochrones from Spain towards the Balearic Islands, thus 
suggesting that late migrants may embark in a more direct cross of 
Mediterranean than early migrants (Figure 3A). 

The maps of autumn migration were also based on a model that 
fitted the data with great accuracy (R 2 = 0.97). They indicate that 
the first 1 5 % of Barn Swallows has already crossed Gibraltar on 1 8 
August (230; Figure 3B). A large migration divide seems to occur 
in France, with Barn Swallows moving along two main migration 
routes, one along the Atlantic coast toward Spain and Gibraltar, 
and the other across Switzerland and along the Italian peninsula, 
with a possible crossing of the Mediterranean from central Italy 
toward Tunisia, thus embarking on a direct Mediterranean 
crossing [33]. A map of the last 15% of migrants also suggests 
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Figure 1 . Progression of Barn Swallow migration in the British Isles. Contour plots of the calendar date in which the CAR model predicts that 
a given percentage of Barn Swallows have been recorded during (A) spring and (B) autumn migration. Contours were generated by linear kriging 
interpolation. Numbers in the colour scale represent the mean date for each 4-days (spring) or 2-days (autumn) colour belt (1 January = 1). For ease 
of interpretation we here report some reference dates: 100 = 31 March, 120 = 30 April; 150 = 30 May, 180 = 29 June, 200 = 19 July, 230 = 18 August, 
260= 17 September. 
doi:1 0.1 371 /journal.pone.01 02440.g001 



that Barn Swallows from north-eastern Europe may move 
westwards across the Balkans and reach central and southern 
Italy, although paucity of data from Eastern Europe prevented a 
clear assessment of movement patterns in this area. 

Phenology estimated by our model was generally consistent with 
that observed in different areas of Europe, as indicated by the 
significant repeatability among observed arrival/departure dates 
and those estimated by our model (Table SI, R = 0.43 ±0.1 6 SE, 
F 2 8 i 97 = 2.51, P = 0.009). Repeatability analyses conduced on 
spring and autumn data separately showed a significant repeat- 
ability between observed and estimated phenology for spring 
migration (R = 0.55±0.16 SE, F 18>19 = 3.49, P = 0.005), and a 
non-significant repeatability for autumn migration (R= 0.25 ±0.33 
SE, F 8>9 = 1.68, P = 0.227). 

Temporal variation in migration phenology 

CAR mixed models restricted to the British Isles did not show 
any significant variations in median date of spring migration, 
either according to period (Likelihood Ratio Test: x 2 2 = 3.19, 



tion (% 2=1-36, P = 0.507, details not shown). Median autumn 
migration date in the British Isles did not change significantly 
between periods (x 2 2 = 0.52, P= 0.770), latitudinal belt {% 2 i = 0.26, 
P = 0.607) or their interaction (% 2 2 = 1.06, P = 0.590, details not 
shown). 

Median spring migration date changed significantly between 
latitudinal belts in western Europe and north Africa (% 2 2 - 21.55, 
P<0.001), but not between periods (x 2 2 = 2.07, P = 0.355). The 
belt by period interaction was non-significant (% 2 i =3A3, 
P= 0.488; Figure 4A). 

The model fitted to autumn data indicated that the belt by 
period interaction was significant (/ 2 4 = 11.81, P = 0.019), as well 
as the main effect of belt {% 2 2 = 20.34, P<0.001). Post-hoc tests 
showed that migration was later in southern than in central and 
northern Europe (z^ — 3.21, P^0.004) and that in southern 
Europe autumn migration post-1990 was 1 3.36± 3.80 SE days 
earlier than pre- 1970. There were no significant differences in the 
timing of autumn migration in central and northern Europe 
(Figure 4B). 



P = 0.203), latitudinal belt {% 4 



.61, P = 0.205), or their interac- 
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Figure 2. Maps of Barn Swallow movements. Each line connects the ring and recovery position of individual Barn Swallows in (A) March-June 
and (B) August-October in the British Isles or in (C) February-June and (D) August-November in western Europe and north Africa. To facilitate the 
interpretation of the figure only Barn Swallows that moved more than 1 and less than 8 degrees latitude or longitude are shown. 
doi:1 0.1 371 /journal.pone.01 02440.g002 



Discussion 

In this study we analysed ring recoveries spanning from 1 908 to 
2008 to describe patterns of bird migration and their long-term 
temporal trends. A first set of analyses was based on a subset of 
ring recoveries in the British Isles, and a second set tentatively 
extended the analyses to Western Europe and North Africa, where 
data are sparser. 

We produced maps describing both spring and autumn 
migration phenology over the British Isles and, tentatively, over 
western Europe and north Africa, from which the main migration 



flyways could be inferred. We found no significant changes in 
migration phenology of the Barn Swallows in the British Isles, but 
an earlier timing of autumn migration in southern Europe and 
north Africa in 1991-2008 compared to 1908-1969. No change in 
autumn migration phenology was observed in central and 
northern Europe as well as in spring migration phenology over 
the whole western Europe and north Africa. 

Ring recoveries currently represent the largest and only long- 
term datasets on bird migration, yet they are hampered by several 
potential sources of bias, primarily due to large spatial and 
temporal variation in sampling effort. Previous studies that faced 
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Figure 3. Progression of Barn Swallow migration in western Europe and north Africa. Contour plots of the date in which the CAR model 
predicts that a given percentage of Barn Swallows have been recorded during (A) spring and (B) autumn migration. Contours were generated by 
linear kriging interpolation. Numbers in the colour scale represent the mean date for each 10-days colour belt (1 January = 1). For ease of 
interpretation we here report some reference dates: 100 = 31 March, 120 = 30 April; 150 = 30 May, 180 = 29 June, 200=19 July, 230 = 18 August, 
260=17 September, 300 = 27 October. 
doi:1 0.1 371 /journal.pone.01 02440.g003 



the same problem tried to account for these potential sources of 
bias, by re-running the analyses on different subsets of data or by 
including additional variables (see e.g. [18;19]). These approaches 
could not be applied in this case, because this study focuses on 
migration periods only, and therefore analyses cannot be restricted 
to 'focal' periods of migration without losing important informa- 
tion on early and late migrants, which are relevant for modelling 
the progression of migration in a given area correctly. In addition, 
results from analyses accounting for other potential sampling 
biases may be difficult to interpret. For example, recovery 
condition of individuals, indicating, for instance, whether an 



individual was actively trapped or fortuitously recovered, and 
whether a bird was dead or alive at the time of recovery, were 
included in previous analyses of ringing recoveries [18]. However, 
trapped birds may provide early- or late-biased estimates of the 
timing of bird migration depending on the scheduled activities of 
ringing stations, which may vary between years and geographical 
regions. Recoveries of dead or live birds may also show different 
biases. For example, if birds are more likely to be found dead early 
than late in spring, and late than early in autumn, analyses 
restricted to birds found alive or dead may bias the outcome in 
opposite directions. In addition, analyses restricted to subsets of 
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Figure 4. Boxplot of median A) spring and B) autumn migration 
dates in western Europe and north Africa. Dates were estimated 
at each cell (t 0 parameter of cloglog curves interpolated at each cell) in 
all belt-by-period combinations. The solid line represent the median 
value, the top and the bottom of the boxes represent the first and the 
third quartile while whiskers approximately include 95% of data. Circles 
represent outliers. Numbers represent sample size (i.e. number of cells 
per period and belt). Asterisk denotes the belt that differs significantly 
from the others atTukey post-hoc tests (z£-3.212, PsO.004 in all cases). 
Different letters denote periods that differ significantly to each other 
within each latitudinal belt at Tuckey post-hoc tests (z = 3.612, 
P = 0.001). 

doi:1 0.1 371 /journal.pone.01 02440.g004 

data may be more prone to produce biased results due to lower 
power of analyses based on reduced sample size than the whole 
dataset, which will give the most robust and the least biased results. 
All the analyses presented in this study therefore used the whole 
dataset, and all the results should be considered by keeping in 
mind that it was impossible to conduct additional analyses 
accounting for possible sources of bias. 

Phenological estimates from our models were generally consis- 
tent with phenology of the spring and autumn migration described 
by both quantitative and qualitative observations in literature. We 
acknowledge that on this analysis we were forced by the scarcity of 



the data to cumulate information from spring and autumn 
migration as well as from quantitative and qualitative observations. 
However, repeatability of observed (from direct observation or 
trapping of live birds) and estimated (from our models) values was 
significant, thus confirming the general consistency of our 
estimates and the observed phenology. Closer observation of the 
data reported in Table SI, however, shows some differences 
between observed and estimated phenology. Comparison with 
time series of first arrival dates indicates that our model estimated 
that the first 5% of Barn Swallows arrive 1 (Norfolk, England [34]) 
to 31 days later (Leicestershire, England [35]) than the first 
swallow was observed. Arrival time of the first 5% of Barn 
Swallows should be close to but later than that of first observations. 
Time shifts of 21-31 days between our estimates and observation 
at sites like Parchim (Northern Germany [36]), Brescia (Northern 
Italy, [37]), or in Leicestershire (Table SI) are therefore not 
negligible, even if we consider that the first Barn Swallows usually 
arrive much earlier than the bulk of migration [38]. 

Estimates of median arrival dates from our model were 
consistent with those of the only published time series of mean/ 
median arrival date that was available to us (V entotene, Southern 
Italy [39]; Table SI). In addition, we were able to reconstruct 
arrival dates of 15%, 50% and 85% of Barn Swallows in Kraghede 
(Denmark) [40] . In this case our model estimated phenology about 
one month later than that observed (Table SI), suggesting that our 
results may depict a somewhat later phenology than site-specific 
time series, at least at some geographical areas. This later estimate 
of timing of migration provided by our results may therefore 
suggest that early migrants are underrepresented in ring-recover- 
ies, that our method underestimates the proportion of early 
migrants, or both. In addition, spring migration spans several 
months and some records may refer to Barn Swallows captured at 
breeding sites some times after that they have arrived. At the same 
time, it is questionable whether time series of arrival at single, 
selected localities can reliably reflect arrival dates at areas as large 
as four degrees latitude per longitude. This latter interpretation is 
suggested by the consideration that time series of arrival dates are 
usually collected at ringing stations, which are located at key places 
along migration routes, and at localities that may be close to the 
margins of a cell. However, our estimates are consistent with the 
general, qualitative description of migration phenology over larger 
geographical areas (Table SI). 

Comparisons of results about the timing of autumn migration 
are more problematic, since information of autumn migration 
phenology is sparser than that on spring migration. Published 
information from Spain [41] depicts an earlier departure of Barn 
Swallows than that estimated by our model. However, our model 
estimated that the last 5% of Barn Swallows are still in central 
Spain (cell G2, Figure SID) on 25 September (268), a date very 
close to that of 21 September reported by Gordo & Sanz [41] as 
the mean departure date of the last Barn Swallow from Spain. 
Mean autumn passage date of Barn Swallows at the Col de 
Bretolet (Switzerland) is 19 days later than median passage date 
estimated by our model (Table SI). Similarly, departure date of 
the last Barn Swallow at four ornithological observatories in the 
UK is 30 to 54 days later than the date at which only 5% of Barn 
Swallows are still in the area estimated by our models. Conversely, 
mean departure dates of Barn Swallows from northern Italy based 
on a short (15 years) unpublished time series (R. Ambrosini, 
unpublished data) are earlier than the date estimated by our model 
for the presence of the last 5 % of Barn Swallows in the cell (Table 
SI). In conclusion, evidence of the ability of our model to 
accurately estimate autumn migration phenology is not unequiv- 
ocal, but information on the timing of departure of 'extreme' 
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individuals at selected localities may not properly represent the 
general timing of migration over larger geographical areas. 

Different independent sources of information however con- 
firmed the reliability of our modelling approach. The pattern of 
spring migration depicted in the maps of the British Isles we 
produced is qualitatively consistent with the south-west to north- 
east pattern of Barn Swallow migration progression through 
Britain described by Huin & Sparks [42] by using phenological 
records compiled before 1947. The fact that our maps indicate a 
migration phenology about 20 days later than that described by 
Huin & Sparks (e.g. 124 to 128 = 6 to 8 May in central England in 
our map vs. 109=19 April in Huin & Sparks [42]) can be 
explained by considering that the maps by Huin & Sparks are 
based on first observation dates while ours are based on the 15% 
of migration movements. In addition both in our maps and in the 
paper by Huin & Sparks the time-difference between isochrones 
though Britain is 20 days (Figure 1A and [42]). Hence, despite the 
difference in the timing of migration, both studies consistendy 
indicated that Barn Swallows arrive in northern Scotland about 20 
days later than in southern England (approximately 800 km to the 
south). 

The second main aim of the present study was to investigate 
variation over time in timing of migration. No significant change 
in spring and autumn phenology in the British Isles was detected, 
which is consistent with Mason [43], who did not find a long-term 
trend in arrival dates of Barn Swallows between 1942 and 1991, 
and with Sparks & Carey [34], who found only a slight trend 
toward a later arrival over two centuries. 

No significant change in timing of spring migration appeared 
neither in the analyses on western Europe and north Africa. Barn 
Swallows are known to have advanced first arrival dates 
throughout Europe. For example, Barn Swallows advanced first 
arrival date by 13 days in 1970-2004 in the Iberian peninsula 
[41], and mean/median arrival dates advanced by 0.34 days 
year" 1 in 1982-2006 in northern Italy [37] and by 0.17 days 
year -1 in 1960-2006 in Europe [8]. Our model was therefore 
unable to capture this widespread advancement, probably because 
the paucity of data from ring recoveries forced us to calculate 
median arrival dates over periods as long as 20 years or more and 
because of a large heterogeneity among periods in the variance of 
median arrival dates in cells (despite accounting for this problem in 
the statistical analyses; see Methods). However, arrival dates of 
migrant birds, and of Barn Swallows in particular, may have 
varied non-linearly over the study period [41]. Indeed, Barn 
Swallows in Spain delayed their arrival dates during the seventies 
and then have advanced, reaching the same arrival dates as pre- 
1970 only in recent years. Our analyses are pardy consistent with 
this pattern. Median spring arrival dates seem to have been 
delayed in 1970-1990 in southern Europe and have then 
advanced, returning to pre- 1970 levels in the last decades 
(Figure 3A). In addition, cloglog curves interpolating arrival dates 
in central Spain (cell G2 of Figure SI) pre- 1970 and post- 1990 
almost overlap, while that for 1970-1990 was shifted towards later 
arrivals (Figure S3). 

Parameters of the cloglog curves indicate a significant advance- 
ment in autumn migration in southern Europe and north Africa 
(Figure 4B). This pattern is consistent with that found by Gordo & 
Sanz [41], who documented an advancement in autumn 
migration in the Iberian Peninsula. In addition, Jenni and Kery 
[44] reported a delay in mean autumn migration at the Col de 
Bretolet in 1970-1982 with respect to earlier years, and a 
subsequent advancement. Also this pattern is consistent with the 
results of our model, that estimated a delay of 6 days in median 
migration date at the cell including this Swiss locality (E5, see 



Figure SI) between 1970-1990 and pre-1970, and a subsequent 
advancement of 1 1 days between 1970-1990 and post-1990 (other 
details not shown). 

In summary, in this study we propose a novel method to 
describe patterns of migrations and main routes followed by 
migratory birds based on ring recoveries. Importantly, this method 
does not use information on the movements of individuals between 
locations where they were observed at different times of their lives, 
but is entirely based on the information on the date at which a bird 
has been observed in a given place. It may therefore be possible to 
extend its application to other, potentially larger, datasets. For 
example, ringing data, which are by far much more abundant than 
ring recoveries, can be used for this purpose. The main 
disadvantage is that - so far - only a few of these data are stored 
in the EDB. They are therefore more difficult to access for 
continent-wide analyses, and they are more prone to temporal 
sampling biases (e.g. non-random variation in sampling effort both 
within and between seasons; [15]). However, they may allow 
detailed studies at the scale of smaller geographical areas (e.g. 
countries). Similarly, this method may be applied to sighting 
databases, such as those collected via the web (e.g. BTO 
BirdTrack project http://bkl.bto.org/birdtrack; the ORNITHO 
family portals e.g. www.ornitho.ch), which are becoming increas- 
ingly popular in recent years, and to databases of timing of 
flowering and leafing [45] . 

Ring recoveries and museum specimens provide the only 
available data spanning over long time periods, and thus they are 
the only data allowing investigation of the variation over time of 
migration phenology over large geographical areas. If ringing data 
too were available over long time periods, the increased amount of 
data available for the analyses may allow the use of reduced 
intervals so that more detailed variation in migration phenology 
over time can be explored. 

Supporting Information 

Figure SI Maps of ring recoveries of Barn Swallows. 

Exact position of ring recoveries of Barn Swallows in (A) March- 
June and (B) August-October in the British Isles and in (C) 
February-June and (D) August-November in western Europe and 
north Africa. 1908-1969: green; 1970-1990: blue, 1991-2008: 
red. Cell ID is shown. Parallels separating latitudinal belts are 
shown. 
(TIFF) 

Figure S2 Example of data interpolation. Complementary 
log-log curve was interpolated to the March-June data from cell 
El in the British Isles (see Figure SI). Dots represent cumulated 
proportion of Barn Swallows recovered in this cell at different 
dates. All data in 1908-2008 were used. The dashed line 
represents June 30. 
(TIFF) 

Figure S3 Interpolated curves showing shifts in migra- 
tion phenology among periods. Complementary log-log 
curves were interpolated to cumulated proportions of Barn 
Swallows recovered during (A) February-June and (B) August- 
November during three different periods at cell G2 in western 
Europe and north Africa (Spain, see Figure SI). Green: curve 
fitted to data in 1908-1969; blue: curve fitted to data in 1970- 
1990; red: curve fitted to data in 1991-2008. In A) the delay in 
spring migration timing in 1970-1990 is evident as well as 
similarity in spring migration phenology pre-1970 and post-1990. 
Curves in (B) evidence the advancement in autumn migration 
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phenology post-1970. The dashed line in (A) represents June 30 

while that in (B) represents August 1 . 

(TIFF) 

Table SI Comparison of observed and estimated 
phenology. General description of timing of migration and 
quantitative information from time-series of arrival dates of Barn 
Swallows in the British Isles and in western Europe and north 
Africa were collected from the literature, websites or other 
unpublished datasets, and compared with the corresponding 
estimate from our models. 
(PDF) 

Text SI The complementary log-log and logistic func- 
tions. A brief description of complementary log-log and logistic 
functions. 
(PDF) 
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